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("-function methods are used to study the properties of the non-relativistic interacting Bose gas 
at finite temperature and density. Results for the ground state energy and pressure are obtained at 
both zero and finite temperature. The method used does not restrict the form of the interaction, 
which can be completely general. A similar procedure is then applied to evaluate the ground state 
energy of a binary mixture consisting of different bosons. Analytical results are obtained in the case 
where the two species have the same mass. 
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I. INTRODUCTION 

The non-interacting Bose gas, and its behaviour as a 
function of temperature, is a key part of any study of 
statistical mechanics. (See for example, Ref. |[ [| 
However the neglect of interactions among the particles, 
which may be justified as a first approximation, cannot 
be valid at any fundamental level. (As emphasised in 
Ref. H , the condensation of a free Bose gas has unphys- 
ical properties, such as infinite compressibility.) Accord- 
ingly the study of interparticle interactions is of obvious 
physical interest. 

The interacting Bose gas received a great deal of at- 
tention about forty years ago. Textbook treatments may 
now be found. (See for example Ref. || [|, J5[.) The 
methods used to study the interacting Bose gas go be- 
yond what would normally be described as elementary. 
The pioneering studies include || [?], |j| who obtained the 
ground state energy of the interacting Bose gas at zero 
temperature using a variety of methods. The generaliza- 
tion to T 7^ was performed by Lee and Yang || . Much 
more recently it was shown how modern functional inte- 
gral methods of quantum field theory could be used to 
study the interacting Bose gas to great effect |lCj . It was 
also shown how the results for the non-relativistic 
interacting Bose gas could be obtained as a limiting case 
of the relativistic gas. 

The main purpose of the present paper is to show 
how it is possible to study the interacting Bose gas at 
low temperatures and densities in a relatively straight- 
forward way. The method makes use of the background 
field method to obtain the excitation spectrum of 
the gas in a manner akin to that of Bernstein and Do- 
delson jl3|]. There is no need to choose a specific form 
for the interaction. The zero temperature result for the 
ground state energy can then be understood as a regu- 
larized zero-point energy and studied in a similar way to 
the Casimir effect |14 . By adopting (-function regular- 
ization IT5L and specialising the form of the interac- 



tion, is easy to derive the classic result || fj], || for the 
ground state energy. Such a technique renders the dis- 
cussion of renormalization of the theory trivial. By an 
examination of the thermodynamic potential it is possi- 
ble to obtain finite temperature corrections to the zero 
temperature results in a straightforward way. 

To further exemplify the methods used we will also 
study a mixture of two interacting Bose gases. Although 
there has been some work on the subject 17, I8[ |l9f so far 
as we are aware, the analogue of the ground state energy 
calculations of Ref. ||, 0, have not been given before. 
We do this calculation for T — as well as T ^ below. 
It is worth pointing out that mixtures of two Bose gases 
have been achieved experimentally for trapped atomic 
gases |20|. (A small sample of some of the theoretical 
work on trapped mixtures includes plL E3, E3l f24j.) 



II. SINGLE SPECIES GAS 

A. General approach 

In order to see how the method of (-function regu- 
larization works in a simple setting we will first study 
the interacting single species gas. In the case where the 
interaction is quartic in the fields it will be shown how 
the classic low temperature results of Refs. 0, 0, ||, ^ 
are found in a relatively easy manner. After describing 
the basic setup in this section, we will obtain th e gro und 
stat e energy and the equation of state in Sees. II B and 
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We begin with the action functional for a single com- 
plex scalar field ^> 

S[*f,^} = J dt J do- B {^(*t*-*t*) 

-J-|W| 2 + M |*| 2 - W 2 )}. (1) 
2m > 

Here £ denotes the space of interest. Our concern here 
will be with flat 3-dimensional space, so that £ can be 
thought of as a finite dimensional box with the infinite 
volume limit taken. The chemical potential /i enforces 
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the conserved particle number. J7(|vl/| 2 ) is the interaction 
potential which we allow to be an arbitrary function. As 
the notation suggests, we assume that U depends on W 
and ^ only through | \£> | 2 , and we take U to be real. 

The thermodynamic potential f2 at lowest order in 
quantum field theory is comprised of three pieces : 

f2 f^class "i" ^T— ^T^O ' (^) 

Here 

^class = ^([/(|§| 2 )-/i|*| 2 ) (3) 

is the contribution from the classical part of the theory 
with W representing any non-zero condensate. (\P can be 
thought of as the background field in the background field 
method We take ^ to be constant here. V is the 

volume of E. The terms in arc the first 

order quantum corrections to the classical theory, with 
^t/o representing the finite temperature corrections to 

the zero temperature theory described by ^yl . The last 
two terms of (||) can be obtained in a consistent devel- 
opment using the background-field method; however, we 
will use a more elementary approach here. 

We will use the method of excitation energies as de- 
scribed in Ref. Begin with the equation of motion 
for 4" and ^ which follow from ([!]) as 

= W + — V 2 ^ + ^ -U'(\t>\ 2 )^ , (4) 
2m 

= + — V 2 * 1 + /if 1 - U'(\y\ 2 )^ . (5) 

2m 

If we expand \P = ^ + $, where $ represents the fluctu- 
ation about the background field then $ satisfies 

= + — V 2 $ + /i$-C/'(|*| 2 )$ 
2m 

-C/"(|*| 2 ) (|*| 2 $ + * 2 $ t ) , (6) 

= + -^V 2 $ t + - C/'(|*| 2 )$ t 

2m 

-C/"(|*| 2 ) (|*| 2 $t + (*t) 2 $) . (7) 

If we define {/^(x)} to be a complete orthonormal set of 
solutions to 

- V 2 / fe (x) =C7 fe / fe (x) , (8) 

we can expand 

*(t,x) = ^Afce-^VfcW , (9) 

*t(t,x) = J2 B ^ lEkt fk(*) , (10) 

fc 

for independent expansion coefficients A k and B k . E k 
are the excitation energies. Because of the assumption of 



completeness for the /fc(x) we find, after substitution of 
(|) and (0) into (|) and (@), 

= ( J E fe --^ + ^-C/'(|*| 2 )_|*| 2 C/"(|*| 2 )) J 4 fc 

V 2m / 

-* 2 {/"(|*| 2 )B fe , (11) 

o = (-E k -^+v-u\m 2 )-m 2 u"m 2 ))B k 

V 2m / 

-(*t)2cr"(|* |2)^ fc . (12) 

For a non-trivial solution for the expansion coefficients 
Ak and B k we must have 

^ = £- M + C/ ' (l * |2) ) V2 (13) 
x(|^-M + ^(l*| 2 ) + 2|*| 2 C/"(|*| 2 )) 1/2 . 

We can now write down the last two contributions to 
the thermodynamic potential ([!]). We have 

n£Lo = 5E*k. (14) 

fc 

since at zero temperature the only contribution is the 
sum over zero-point energies. Also 

^T 1 io = ^E ln (l-e- /5£fe ), (15) 

k 

which is the expression usually written down for the ther- 
modynamic potential at finite temperature. We leave the 
factor of h in front of both expressions ( fbi]) and (|l5[) to 
denote that they are of first order in a systematic loop 
expansion. We will only work to order h. 

^t=o * s divergent, and to make sense of this we must 
regularize it in some way. We will use a version of £- 
function regularization |L5], |l6| here. (See Ref. J2^] for 
a pedagogical introduction to the method we will use.) 
This consists of defining an energy (^-function by 

E{ S ) = Y. Ek ( lEk )~ s > ( 16 ) 

k 

where s is a complex variable suitably chosen so that the 
sum in ( |l6| ) converges, and £ is a constant with dimen- 
sions of inverse energy introduced so that E(s) has units 
of energy for any value of s. The regularization consists 
of defining Clj^n by 

4=o = ^(0) , (17) 

where E(0) represents the analytic continuation of E(s) 
from the region of the complex s-plane where (]l^) con- 
verges to s = 0, assuming that E(s) is regular at s = 0. 
It may be that E(s) has a simple pole at s = in which 
case E(0) in ( |l7| ) represents the analytic continuation of 
E(s) to a neighbourhood of s = 0. The pole can be dealt 
with by the usual machinery of renormalization. In the 
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theory considered here renormalization is not necessary 
as we will show that E(s) is analytic at s = 0. 

Given the thermodynamic potential Q, the number of 
particles is given by 



N = - 



and the internal energy density p is 



Using (|18[) it is easily seen that 



In addition, the condensate ^ satisfies 

dfl 



= 



(18) 



(19) 



(20) 



(21) 



The exact equation for ^ is very complicated due to the 
dependence of the excitation energies E k on 'J in (13). 
However because we are only working to order h, from 
(0) and (H) we can conclude that 



p = u\m 2 )+o{h) 



(22) 



This means that in any expression which is already of 
order h we can simplify it using only the first term of 
( p2| ) . Of course in any term which is of lower order in h 
we must include the 0(K) corr ectio n to \i for consistency. 
We will encounter this in Sec. II C. 



B. Ground state energy 

We will introduce the number density n by 

N 

n = — . 
V 



(23) 



From (||,| 15 T^) we find the number density of particles 
to be 



n = |*| 2 +4!1 + n^ 



with 



l T=0 



1_ ( dfl ( ^l " 
V [ dp 



(24) 



(25) 



and a similar expression holding for rSpL Q . Since rSp =0 



and of order H, it is obvious from ( |24| ) that 



|*| 2 = n + 0(H) . 



(26) 



Using the same argument as for p below (|22|), we can set 
l^l 2 = n in any expression which is already of order %. 

If we use (|2]j|) along with (|2^) and (|i|) in the expres- 
sion (EG) for the energy density, we have 



1 JL 

V 



(27) 



We can use (|2J) to simplify the argument of the energy 
J7(|>I , | 2 ). To order h we have 



C/(|§| 2 ) = U(n) U'(n)(4l +n^>) + 



(28) 



The last term of (| 
(p6|) to give 



may be simplified by using ( |22|) and 



K n T ] =o + n %o) = u '( n )( n TL + 4yo) + • • • > 
again working only to order h. We therefore find 



(1) 



U{n) 



1 

y~ 



-n 



(i) 

T=0 



1A 

ya/3 



(29) 



(30) 



from fl27|), holding to order ft. If we use (|15|), then the 
last term of ( |30| ) becomes 

1 d 



o {1) 



E k (e 



iy 



(31) 
(32) 



This is the usual expression for the internal energy den- 
sity from statistical mechanics. 

We now require the evaluation of £1^=0- ( P^l ) can 
be used for this. A direct evaluation using the energy 
levels in (|l3|) is difficult. Because we are only working to 
order h, and ^^1 is already of order h, we can simplify 
the excitation energies using p — U'(n) and l^l 2 = n. 
This results in 



E k 



\2m) V2m 



2nU"(n) 



1/2 



(33) 



If we take the large volume limit we can replace Gk — * k 2 
and Ylk ~ * ^ / (2^)3 i thus the energy (- function to be 
evaluated is 



V 



d 3 k 



*1 
2m 



(l-s) 



— + 2nf7"(n) 
2m 



(34) 

(Note that the energy spectrum in Eq. ( |33| ) with o& — > 
fc 2 was first given by Bogoliubov 26 .) The integral in 
(^) is easily evaluated using the representation of the 
T-function to give 



E(s) = (2m) 3 / 2 



y r(2-f)r( s -f) 



r(^) 



[2n?7"(n)] 5 / 2 - 



(35) 
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This expression is analytic at s — with the result 
16 



E(0) 



15tt 2 

We then find from ( |l7| ) that 



m 3 ' 2 V[nU"{n)fl 2 . 



(36) 



li Lrp n 



8/i 



V T=0 ~ 15^ 2 
If we use (B2l) we have 



m 3 ' 2 [nU"{n)f/ 2 



d 3 k 



El; 



(2tt) 3 (ef 3E » - 1) ' 



(37) 



(38) 



where is given by (|33|) after the large volume limit has 
been taken. Using spherical polar coordinates, perform- 
ing the angular integration, and rescaling k — > (2mT) 1//2 fc 
results in 



o {1) 



2n 2 (3 \ (3 



2m 



3/2 



dkk 3 ^k 2 + 2[3nU"{n) 



,k^k 2 +2/3nU"(n) _ ^ 



(39) 



Although it is not possible to evaluate the integral in 
( |39| ) exactly, it is relatively straightforward to find the 
low (or high) temperature expansion. We concentrate 
on the low temperature expansion here and assume that 
2(3nU"(n) >> 1. The first two terms in the asymptotic 
expansion of (39) give 



(1) far 2 / to \ 4 r 
Pt ^° " 30 J ^ 

+ ...}. 



25tt 2 
21 



[/?n[/"(n)r 



(40) 



We then have the low temperature expansion of the en- 
ergy density in (BOh being 



p(T) ~ (7(n) + 
far 2 T 4 ' 



30 



8fi 
15tt 2 
m 



m 3 / 2 [n(7"(n)] 5 / 2 



(41) 



nU"(n) 



3/2 



25tt 2 
~2T 



[/3nt/"(n)]" 



In the zero temperature limit only the first two terms 
survive and we have the exact result 



p(T = 0) = U(n) + 



8h 

15^ 



m 3 / 2 [nU"{n)f/ 2 



(42) 



So far we have kept the potential C/(|*| 2 ) arbitrary. In 
order to make contact with the standard results we will 
take 



um 2 ) = 



2ixa , 



1*1 



(43) 



with a the s-wave scattering length. In this case it follows 
from (E§) that 



p(T = 0) 



2-Kan 1 



m 



1 



1287i 



(44) 



and from (|40|) that 
hm 3 T 4 



Pt=o 



240^ 



(na 3 ) 



3^,-3/2 



! 25 /toTV 
336 \ na J 



(45) 

The result in (44|) agrees with the classic calculations 
of Refs. P, 0, |§ for the zero temperature gas. The first 
term of (E9) reproduces the leading correction at low tem- 
perature found in Refs. ||. (Note that Ref. || chooses 
to = 1/2.) 



C. Equation of state 

The equation of state involves a determination of the 
pressure, which follows from the first law of thermody- 
namics as 



P = - 



dV 



(46) 



Because Q is an extensive quantity, in the large volume 
limit we have simply 

P=-y- (47) 

The evaluation of the equation of state involves the de- 
termination of Q to order fa From (Q) and (|^) we have 



q(1) i 

T— T^O 



(i) 



1*1 



u(\m 



V 



(48) 



Using (j24|) and 
h) 



1) to eliminate 1*1 2 we find (to order 



P = pn — U(n) 



(i) 

T=0 



(1) 
T^O 



V 



(49) 



where we can set | * | = n in since these 

two terms are already of order fa What makes the calcu- 
lation of P more difficult than that of p is that we cannot 
simply use ( [22] ) for p in the first term; we need to know 
the O(h) contribution to p coming from the quantum 
corrections to the classical field equation for 

Because we are interested only in the condensed phase, 
we can replace (^lj) with 

o=' dn 



This leads to 



p = u'm 2 ) 



9|*| 



1 0(4=o 



P,V,fi 



+ 4io) 



V 



9|*| 



(50) 



(51) 



We now use ( pjj ) again in the first term to find 

i / 0(4=o + 4%) 



V 



+U"{n) 



am 



(1) 

T=0 



(i) ^ 



(52) 
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Using this in (Jf9|) results in 

1 f 9ft (1) 
P = nU'(n)-U(n) + v ln-^ 



+nU (n)- 



+nU"(n)- 



(i) 

T=0 



1 

V 



( =-i 

J d|*| 2 



7V0 



djjL 



(i) 

T#0 



(53) 



for the pressure to order h. The last term of (|53j) displays 
the finite temperature correction to the zero temperature 
result. In the last two terms of ([53|) we can set |*| 2 = 
n and fi = U'(n) to simplify the results after we have 
performed the indicated differentiations. 

We regularize the terms in £1^=0 by the C~f unc ti° n 
method as described in Sec. II A. From the excitation 
energy (Oh it is easy to see that 



dE k 1 / a k \ 
— = I h nu (n) J 



E k V2m 



(54) 



Hp = i{( 2t/ » + ^») 

+n(U"{n)) 2 } . 



2m 



(55) 



It then follows that 



V I 91*1 



+ U"{n) 



on 



(i) ' 

T=0 



x [U"(n) + nU'"(n)}J2E,: 1 - s ^{56) 

k 



with the right hand side analytically continued to s = 0. 
Taking the large volume limit it is easy to show that 



n I 9f2 



(i) 

T=0 



V 9|*| 



+ U"(n) 



on 



(i) ' 

T=0 



hn 



dji 



(2m) 



3/2 



x[U"(n)+nU"'{7i)]r s {l-s) 

r(2-*/2)r(*- 3/2) ,„_„,„_ 3/2 _ 5 



-[2m£/"(n)r 



= -^[[/» + n(/"'(n)][mnC/"(u)] 3 / 2 (57) 

after analytic continuation to s = 0. 

The pressure of the zero temperature gas follows from 

P(T = 0) = nU'(n) -U(n) + ^m 3 / 2 \nU"(n)} 5 / 2 

+ ^-n 2 U'"(n)[mnU"(n)] 3 / 2 (58) 
3tt z 



if we use ( p7| ) and <pi[). Specialising to f/(n) = 2-Kar? jm 
as in (^H) results in 

27ran 2 



-Pt=o = 



64/i 



(59) 



in agreement with the standard results g, raj- 

To evaluate the finite temperature part of the pressure, 
we use ( p"5| ) along with (Q) and (55) to find 

- £[tf>) + »tf'>)] 



an 



V 9|* 



-1 



V 



1 1/2 



(Tfe + 4mnU"(n) 



(60) 



In the large box limit we can evaluate the right hand side 

at low temperatures as we did for /9yi m ^ ec - A 
short calculation leads to 

v \ 9|*| 2 + ( ' dn 

nU"'(n) 



hn 2 T 4 



■m 



60 [nU"{n) 



3/2 



1+ wJ{ 1 -Vi 3 " f '"<")]- 2+ -j < 61 > 

where (3nU h '(n) >> 1. In a similar way the low temper- 



ature expansion of ( |15[ ) gives 

TT 2 h 



n (1) 



V 



90 



nU"(n) 



3/2 



-^W>)r 2 



} • (62) 



Combining © and (||) we have as the finite tempera- 
ture correction to the pressure 



36 



" rp4 



m 


3/2 | 


nU"(n)_ 





3nU"'(n) 
5U"{n) 



17" (n) 

at low temperatures. When we specialise to U(n) 
2-Kar? j m this result reduces to 



63) 



P 



^m 3 T 4 (n«r 3/2 |l-— f^V + 
288 I H8 \(3na J 



. (64) 

The leading term agrees with the earlier calculation in 
Ref. §. 



III. BINARY MIXTURE 

A. General approach 

We now generalize the results of Sec. II to a two-species 
Bose gas described by Schrodinger fields *i and * 2 . We 
allow for self-interactions as well as mutual interactions 
between the two species in a general way. The action 
functional will be chosen to be 

2 



s = 



*/*.[£{!(•}♦,- 

S 3=1 

^|*/}-t/(|*i| 2 ,|* 2 | 



i 



2m j 



IV*,- 



(65) 
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Here toi j2 denote the masses of the two species, and 
are the chemical potentials. We assume that the inter- 
action, described by the potential U, only involves the 
magnitudes of the two fields. The theory has a rigid 
U(l) x U (1) symmetry in general, unless we further spec- 
ify the potential. The field equations which follow from 
(f35|) are 



= fflj 



_Lv 2 * 7 

2m j 3 



l*2| 2 



(66) 



for j — 1,2 along with the hermitian conjugate. Here 

The first requirement is to identify the excitation ener- 
gies by first setting \&j = ^>j + &j and linearizing in 
As in Sec. II A, tyj denotes the background field which 
we assume to be constant. Wc find 



= i$i 



1 



2mi 



V 2 $i + - (CTi + |*i| 2 i7ii)$i 



-* 2 C/n$ t 1 - *i^[/ 12 $ 2 _ ^ 2 £/ 12 $ 2 , (67) 

along with its hermitian conjugate, and two other equa- 
tions obtained by interchanging the labels 1 and 2. We 
have abbreviated Ujf = d ^ - p U = Uj/j. All terms 

involving the potential are evaluated at a background 
fields. To get the excitation energies we set 



x) = j2 A ^- iEkt f k ^) 

k 



(68) 
(69) 



for independent coefficients Akj and Bkj- (/fc( x ) is still 
the solution to (g).) After some calculation the excitation 
energies can be shown to be 

eL = \(<Z + 4-\*i\*v&-Wvk) 



2 



16|^i| 2 |§ 2 | 2 C/ 2 2 



x(a 1 -|* 1 | 2 [/ 11 )(a 2 -|4' 2 | 2 C/ 22 ) 



1/2 



(70) 



where 



2m 



17,- 



(7 = 1,2). (71) 



The thermodynamic potential has the form (|l|) 
where 

^class = ^[f/(|*l| 2 ,|* 2 | 2 )-Ml|*l| 2 -M2|* 2 | 2 ](72) 



n 



n 



(1) _ 

T=0 ~ 
(1) 



T#0 



ar£[in(i 



-/3-Ej, 



■lnfl 



-£.E*- 



)]• 



(73) 



(74) 



The number density for species j is (see dig)) 

1 /an 



(J = l,2) 



(75) 



T.V,* 



and the energy density is given by (|l9|). It is easy to show 
that in place of ( po|) we have 



1 5 

The pressure is given by (|47j). 

B. Ground state energy 



(76) 



Using (fra) we can write 



n, = I*, 



n 



(i) 



(7=1=2) (77) 



as in (242a). This allows us to conclude that 



Tlj 

we find 

P = 



0{h). Since only fiyi has an explicit j3 dependence, 



cW,|* a | a ) + ^0 + ^(^0) 



+Mi(ni - |*i| 2 ) +/i 2 (n 2 - |* 2 



(78) 



To order h we have fj,j = Uj as in (p2|). We can simplify 
( fr8l ) to order ft as 



U(m,n 2 ) 



(79) 



where 



PT^O 



£7, 



(e^--l) *} (80) 



is the finite temperature correction to the zero- 
temperature result. 

We were not able to obtain analytical results for the 
ground state energy density except in two special cases. 
One case is Un = 0, which means that the two species do 
not interact with each other. This is not very interesting 
since it just results in a sum of terms for the two species, 
with each term given as in Sec. II. The second case where 
we can obtain exact analytic results at T = occurs 
if we take mi = m 2 . A physical system where this is 
relevant concerns different hyperfine spin states of a given 
gas. This specialisation allows a sufficient simplification 
of ( |70| ) to enable analytic results. 

We are free to set | | 2 = rij in the second two terms of 
([79|), as well as to take /ij — Uj. With these two simpli- 
fications made, and mi = m 2 = m taken, the excitation 
energies ([70|j7l]) become 



E k± 



Ok ( Ok 

2m \2m 



a± 



(81) 
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where 



IV. DISCUSSION AND CONCLUSIONS 



a± = niUu + n 2 U 22 

± [(mCTu - n 2 U 22 f + 4nw2Uf 2 ] 1/2 . (82) 

The result for &>t=o ma y be defined using the energy £- 
function method exactly as in Sec. II. In fact the energy 
C-function for the energies E^± may be obtained from 
( p5| ) by replacing 2nU" with a±. The energy density of 
the zero temperature gas is therefore 



Pt=o = U(m,n 2 ) + 



30tt 



■(2m) 



3/2 / 5/2 



5/2 



(83) 



The finite temperature correction, at low tempera- 
tures, can be obtained from ( fiof ) with the replacement 
2nU" with a±. The leading term at low temperature is 



Pt^o 



h—T 4 

30 



2m \ 
a+J 



3/2 



2m \ 



3/2' 



(84) 



If we set U\ 2 = 0, then the energy density reduces to the 
sum of the single species gas contributions. The presence 
of mutual interactions of the two species complicates the 
result. 

It is possible to generaliz e the method described for 
the single species gas in Sec. QC to deal with the binary 
mixture. However as the details are rather messy, and 
the result not particulcarly simple we will not pursue 
this here. 



We have shown how C-function techniques may be 
used to obtain the thermodynamic properties of a single 
species as well as a two species Bose gas. The interac- 
tion in either case was arbitrary. For the binary mixture 
analytical results at zero temperature were obtained in 
the case where the masses of the two constituents were 
equal. Although this restriction might seem artificial at 
first sight, and perhaps not of great physical relevance, 
it has been realised experimentally for trapped gases [p0[. 
Some progress can be made in the more general case of 
unequal masses, but analytical results do not appear pos- 
sible. 



The overwhelming interest at the present time is in the 
behaviour of trapped Bose gases. In the present formal- 
ism this entails the introduction of an additional term to 
(|lf) to represent the trapping potential. This loses the 
feature of the calculation above which enabled analyti- 
cal results to be obtained; namely, that "J/ was constant. 
Nevertheless the background field method works whether 
\& is constant or not, and it seems likely that the use 
of generalised ^-function methods like those presented 
above can be used to provide a different approach to the 
problem of trapped interacting gases. 
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